r.i: ", 


ZA-532 

Fluid  ^  «»iniCB  Research 


0 


f  ^ 


0  O 


""c 


4sr  *a 

"*»C  *«, 

'*■  „*«»<*,  / 
4,  *  * 


ON  THE  MIXINO  FROBUSM  OF  AN  AXI-8WETRIC  FREE  JET 
INTO  AIR  INCIUDING  CHEMICAL  REACTIONS 

Inge  L.  RyhRing 


March  1961 


PROJECT  DEFENDER 


CONVAIR  RAN  DICOO  CONVAIR  DIVISION  OKNKRAL.  DYNAMICS  CORPORATION 


BEST 

AVAILABLE  COPY 


OH  T*  MRUS  1 
HID  AIR 


i  h  n .  i  ii  r 


ah  m-mmaac  mat  jit 

CMOCAL  RkACTIOHB 


Infli  L.  Rytalng 
March  1961 


Tina  raaaarch  «m  aopporttd  Iqr  Admnecd 
Hrojaeta  Agency  Ordtr  Ho.  126-60  (Preset 
Dafaaiar  through  Contract  AF- 19(604) -555V 


COHVAIR 

A  Divlalon  of  Qahcral  fynanlca  Corporation 
San  Ditto 


.■ 


ii-.jr  •. 


Oil  THE  Nixmo  PROBLEM  OF  AX  AXI-8YI§tETRIC  FREE  JET 
mo  AIR  HCUUDXH  CHBGCAL  REACTIOXS 

Inge  L.  Ryhmlng 


ABSTRACT 

y 

Velocity  end  temperature  field*  together  with 
specie  concentrations  are  confuted  for  an  axially 
Bjpuitrl  r,  supersonic,  hot,  free  jet  nixing  with 


quiescent  air  •  Based  on  simplifying  assumptions 
regarding  the  chemistry  involved,  the  problem  can 
be  divided  in  tvo  parts.  One  of  these  parts  is  to 
determine  the  velocity  and  temperature  field.  An-* 

l^The  second  part  con¬ 


sists  of  integrating  the  continuity  of  the  specie 
equations  when  the  solution  to  the  first  part  of 
the  problea  is  obtained.  Two  different  cases  are 
considered  fro  sen  and  equilibrium  flow^.  For 
frozen  flow  it  is  shown  that  the  specie  concentra¬ 


tions  are  similar  to  the  velocity  distribution. 


For  equilibrium  flow  the  specie  concentrations 
are  obtained  with  a  crude  but  slmpla  analysis,  in 
which  the  effects  of  diffusion  and  chemical  reac¬ 


tions  are  combined. 


DURODUCTIOH 

The  gas  dynamics  ,of  mixing  phenomena  in  plane  or  axi-symnetric  free 
jets  into  air  received  considerable  attention  in  the  past.1.  The 
theory  has  been  developed  by  using  the  usual  boundary-layer  approxima¬ 
tions  to  determine  laminar  as  veil  as  turbulent  mixing. 
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If  the  chemical  phenomena  vhich  arise  in  mixing  hot  gaaea  of  vari¬ 
ous  chemical  composition  are  included  in  the  analysis,  the  gas  dynamics 
of  the  Jet  mixing  problem  becomes  much  more  complicated.  The  difficulty 
arises  mathematically  from  an  increase  in  the  number  of  dependent  vari¬ 
ables  and  also  the  development  of  nee  equations  added  to  the  problem, 
e.g.,  the  continuity  of  specie  equations  are  coupled  with  each  other 
and  with  the  energy  equation  of  the  reacting  mixture. 

At  the  present  time  very  little  is  known  about  the  rate  constants 
in  the  production  tents  of  the  specie  equations.  So  f*r,  the  two  limit¬ 
ing  cases  (frozen  flow  or  zero  reaction  rate,  equilibrium  flow  or  infinite 
reaction  rate)  or  small  deviations  from  them  have  in  general  been  con¬ 
sidered.  Also,  the  chemistry  involved  has  usually  been  confined  to  two 
component  mixtures.8 

In  this  investigation  the  velocity,  temperature,  and  specie  concen¬ 
tration  distributions  of  an  aid-symmetric  and  supersonic  multi-component 
hot,  free  jet  jiixing  with  quiescent  air  is  computed  for  frozen  and 
equilibrium  flow.  The  mixing  process  is  supposed  to  occur  under  con¬ 
stant  pressure,  so  that  the  pressure  of  the  Jet  at  the  nozzle  exit  is 
completely  balanced  with  the  surrounding  air  pressure.  In  the  analysis 
several  rather  crude  assumptions  are  made  to  facilitate  a  simple  solu¬ 
tion.  The  simplifications  Imposed  are  mainly  connected  with  the 
chemistry  Involved. 
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Consider  a  gas  which  is  formed  in  the  combustion  chamber  of  a  rocket 
motor  and  which  then  expands  through  the  rocket  nozzle.  Depending  on  the 
mixture  ratio  of  the  rocket  propellant  (RP-1)  and  oxygen  and  on  the  expan¬ 
sion  process  through  the  nozzle,  the  gas  will  have  different  composition, 
pressure,  temperature,  and  velocity  distribution  at  the  nozzle  exit.  Two 
different  mixture  ratios  and  two  different  expansion  processes  will  be 
considered.  First,  it  will  be  assumed  the  expansion  follows  a  frozen 


state  with  a  mixture  ratio  of  6  and  second,  the  expansion  is  followed  by 


an  equilibrium  state  with  a  mixture  ratio  of  2.5*  Furthermore,  the  expan¬ 
sion  is  supposed  to  be  carried  so  far  that  the  pressure  at  the  nozzle  exit 
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is  balanced  completely  with  the  pressure  of  the  surrounding  still  air 
at  sea  level  or  1  atm.  We  also  assume  that  the  velocity,  temperature, 
and  concentration  of  the  various  species  of  the  gas  have  a  constant 
value  over  the  nozzle  exit. 

The  gas  will  be  regarded  as  being  composed  of  the  following  species: 
HeO,  CQa,  CO,  Re,  02,  and  a  small  amount  of  various  other  species  which 
can  be  taken  as  non-reacting.  As  the  temperature  range  we  are  consider¬ 
ing  Is  below  2600 *K  we  do  not  take  Into  account  any  dissociation  of  Qe 
and  Re*  The  reactions  which  may  take  place  are  the  formation  of  CQe 
and  Re  from  CO  and  ReO,  and  the  formation  of  ReO  from  He  and  Qe.  These 
reactions  are  known  to  proceed  according  to  a  number  of  elementary  steps. 
Here,  1  'wever,  the  simplest  possible  form  of  chemical  reactions  will  be 
assumed. 

The  air  outside  the  Jet  is  supposed  to  be  composed  of  Qe  and  H2; 
the  latter  Is  regarded  as  an  inert  gas.  The  temperature  of  the  air  Is 
assigned  the  value  Tfc  *  500 *K. 

A  constant-pressure,  laminar-mixing  process  will  first  be  assumed, 
but  later  on  In  this  investigation,  -turbulent  mixingt*ill  also  be  con- 

S'  - - 

side red.  For  calculating  the  properties  of  the  described  flow  model, 
the  equations  of  motion  of  the  reacting  mixture  written  with  the  usual 
boundary-layer  approximations 'will  suffice.1'3  Thereby,  the  effect  of 
the  temperature  gradient  on  the  speed  of  diffusion  will  be  neglected. 

By  introducing  an  axi- symmetric  coordinate  system  with  the  axis  of  the 
Jet  being  taken  as  the  x-axis,  the  governing  equations*  are: 

continuity  of  species: 


#For  the  exact  equations  of  motion  of  a  reacting  mixture,  see 
Hirschf elder,  et  al.4 
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over-all  continuity: 


JL  (pur)  +  ~  (pur)  =  0, 
dx  or 


momentum: 


energy: 


w  li  w  mm 


Here  p,  u,  v,  and  p  have  their  usual  meaning,  CA  is  the  mass  concentra¬ 
tion  of  species  i,  D1  the  diffusion  coefficient  of  species  i,  and  W  i6 
the  production  function  of  species  i  in  the  chemical  reactions  considered. 
In  the  energy  equation,  h  stands  for  the  enthalpy  of  the  mixture 


h  Ci  ^  ’  hi  ^  Ci  Hl 


where  h/  is  the  heat  evolved  in  the  formation  of  component  i  at  0°K  per 
i  (*_  n 

unit  mass,  and  Pr  is  the  Prandtl  number  (Pr  =  -£-),  and  Lcj  is  the  Levis* 

*  CnODf .  * 


Semenov  number  for  species  i  (Le^  =  ^  ) 


The  system  of  Equations  (1)  through  (4)  is,  however,  far  too  com¬ 
plicated  for  a  simple  solution  to  bring  out  the  main  features  of  the 
mixing  problem.  Further  simplifications,  therefore,  are  needed.  One 
assumption  which  is  generally  employed  and  which  will  be  adopted  here 
is  that  the  Lewis-Semenov  number  equals  one.  Then  the  energy  equation 
is  decoupled  from  the  specie  equations  and  a  solution  to  Equations  (2) 
through  (4)  can  be  obtained  first.  This  simplification  is  justified 
because  Lej^  has  a  value  close  to  one  in  any  actual  case  and  because 
its  variation  for  each  specie  is  small  in  the  flow  field. 


The  equilibrium  flow  case  was  chosen  because  the  chemical  reactions 
which  occur  do  not  play  a  dominant  role  in  the  problem.  Such  a  situation 
can  be  accomplished  by  choosing  a  low  value  for  the  mixture  ratio  of 
propellant  to  oxygen.  This,  in  turn,  implies  that  the  effect  of  reac¬ 
tions  on  the  specie  concentration  will,  in  general,  be  small,  since  the 
largest  part  of  the  reactants  participating  in  the  burning-type  reactions 
already  have  reacted  in  the  combustion  chamber  and  given  non-combustible 
products.  Consequently,  over  the  main  part  of  the  flow  field,  the  W  term 
in  Equation  (l)  will  be  small  compared  to  the  other  terms  involved  except 
possibly  for  a  thin  layer  in  the  mixing  region  where  burning  reactions 
can  be  stimulated  by  oxygen  present  in  the  surrounding  air.  Effectively, 
therefore,  the  equilibrium  flow  field  will  not  differ  significantly  from 
a  corresponding  frozen  flow  field  where  no  reactions  occur  and  the 

W  term  is  zero. 

Since  the  chemistry  is  of  minor  importance  compared  with  the  fluid 
mechanics  of  the  problem,  it  seems  reasonable  to  simplify  the  energy 
equation  further.  Hence,  it  is  assumed  that  the  enthalpy  throughout 

can  be  approximated  as 

h  =  cpT 

vhere  c„  is  .  mean  value  for  the  mixture.  In  this  rather  crude  assump¬ 
tion  ijis  neglected  entirely!  l.a.,  the  contribution  to  the  enthalpy 
caused  by  the  release  of  energy  through  reactions.  Houever,  for  our  equl 
librium  flow,  the  energy  release  will  be  small  because  of  the  chosen 

mixture  ratio. 

By  applying  these  assumptions  to  the  energy  equation  and  the  further 
assumptions  that  the  density  variation  with  temperature  follows  the  ideal 
gas  law  (pressure  constant) 


_e_  _Io  (5) 

Po  "  T 

and  that  the  Prandtl  number  equals  one,  the  problem  is  essentially 
divided  into  two  parts.  The  first  of  these  parts  is  to  find  the 
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velocity  and  temperature  variation  in  the  flow  field.  When  these 
variations  are  determined  the  second  part  consists  in  finding  the 
specie  concentrations  for  the  two  different  cases  we  are  considering. 
The  first  part  of  this  problem  has  been  solved  by  Pai.5  However,  the 
essential  steps  in  the  calculation  are  needed  later  on;  therefore, 
they  are  repeated  here  for  convenience. 


THE  VELOCITY  AND  TEMPERATURE  FIELD 

With  the  assumptions  made  in  the  last  paragraph,  the  energy 
equation  reads 

■  <6> 

As  shown  by  Crocco®  and  others,  a  solution  to  Equation  (6)  is 

T«A+Bu (7) 
^CP 

where  A  and  B  are  constants  determined  by  the  initial  conditions. 

The  following  non-dimensional  quantities  are  now  introduced  in  the 
over-all  continuity  and  momentum  equations: 


x* 


where  index  o  refers  to  the  nozzle  exit.  Then  the  independent  variables 
in  the  momentum  equation  are  changed  to  x*  and  ♦,  where  the  stream  func¬ 


tion  ♦  is  defined  such  that 


=  -  'r*p*v*,  \| =  r*p*u*.  (9) 

The  result  is 


du» 

dx*  ro  Po  u0 


d  (  _  _  _  -2  du 
^*p*u*r*  w). 


(10) 
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The  over-all  continuity  equation  is  automatically  satisfied  by  introduc¬ 
ing  the  stream  function  Equation  (9).  In  Equation  (10)  any  temperature 
dependence  of  p  with  T  can  be  chosen.  Here  this  temperature  dependence 
is  made  as  simple  as  possible  by  taking 


H*p*  =  1. 


Furthermore,  we  introduce 

*o 

i  =  x* - 

roPoUo 

as  a  new  independent  variable.  The  momentum  equation  then  is 


(11) 


(12) 


de  ~  d*  \  ^  d*  /  *  (15) 

For  the  numerical  integration  of  Equation  (1?)  the  behavior  of  the 
momentum  equation  on  the  x-axis  is  needed.  This  condition  can  be  obtained 
from  Equation  (3)  by  taking  the  limit  as  r  approaches  zero  and  by  using 
the  symmetry  conditions 


v(x,o)  *  0,  -r—  u(x,o)  =  0. 

o  r 


The  result  expressed  in  the  non-dimensional  quantities  is 


on  r*  =  0:  u*  =  2T*2  ^u# 


as 


dr*2 


(14) 


In  view  of  the  last  paragraph  the  initial  conditions  for  u*  and  T* 


are 


u*  =  1,  T*  =  1  for  r*  <  1 
Ta 

u*  =  0,  T*  =  m“  for  r*  >  1. 


(15) 
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Equation  (7)  with  Equation  (15)  gives  now,  in  general, 


In  Table  1  temperatures  and  concentrations  of  the  various  species  at  the 
nozzle  exit  are  shown  for  the  two  cases  investigated.  Prom  this  table 
and  Equation  (16)  the  temperature  variation  in  the  flow  field  is  given 
as 

Frozen  flow  T* 

Equilibrium  flow  T* 

Equation  (13)  can  now  be  integrated  numerically  with  the  help  of 
Equations  (l4),  (15),  and  (17)*  In  the  numerical  procedure  it  is  further 
necessary  to  have  the  relation  between  the  i|r  and  r*.  From  Equation  (9) 
this  relation  is 

r*2  =  2  f  d  ♦ . 
jo  u* 

The  problem  was  programed  for  an  IBM  704  computer  and  results  from  the 
computation  cure  shown  in  Figures  1  and  2. 


-  0.1652  +  0.8348  u*  +  U-(;L:,U*J, 

1.0517 


«  0.1161  +  0.8839  u*  + 


THE  CALCULATION  OF  SPECIES  IN  THE  FLOW  FTET.n 

For  calculating  the  various  species  in  the  flow  field  Equation  (1) 
should  be  used.  First  we  transform  it  to  the  x#,  4,  variables  with  the 
result 


dCi* 

dx  * 


■  4  37  (r*2p*2u#D i*  ® + pt 


^  ci. 


(18) 
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FIGURE  2  DISTRIBUTION  OF  VELOCITY  U*  AND  TEMPERATURE  T* 
IN  THE  EQUILIBRIUM  FLOW  CASE 
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Ci  Di 

where  C<*  =  - —  and  D.#  =  - —  . 

1  Ci«  1  Din 


c_p  D. 


But  Le^  =  -*—■ 


=1,  Pr  = 


1  and  c  -  const. 
P 


Hence,  Dj*  p*  =  p#  and  since  (i*  p*  =  1  we  have  p*2  D^*  =  1. 
In  terns  of  Uq,  pQ  and  pQ  we  have 


rouo  ro  Po  ^ 


Thus,  the  same  independent  variables  can  be  used  as  in  the  preceding 
paragraph,  and  Equation  (18)  can  now  be  written  as 


d  /  ,  SC1*\  .  •  ro2 

is  'Si-V1*1*  ■ 


For  frozen  flow  there  are  no  reactions  taking  place  nnd  the  produc¬ 
tion  term  in  Equation  (20)  is  zero.  Equation  (20)  then  has  exactly  the 
same  form  as  Equation  (13)  since,  as  we  cure  using  non-.limensional  quan¬ 
tities,  it  is  immaterial  if  we  speak  about  u*  or  C^*.  The  distribution 
of  a  particular  specie  of  the  jet  is  similar  to  the  u*  distribution, 
because  for  the  jet  species  the  initial  conditions  are 


Cj*  =  1  for  r*  <  1 


Cj*  =  0  for  r*  >  1  . 


We  have  Cj*  =  u*  or  =  Ci  u# 

o 


For  the  0r  and  Nr  concentration  of  the  ambient  air  the  initial  conditions 
are 


V 

=  1 

for 

r* 

IV 

H* 

V 

-  0 

for 

r* 

1 

(«?) 


For  these  species,  therefore,  a  new  dependent  variable 


R 


is  chosen  and  we  have  then  R  =  u»  -  1  -  c  *  or  C.  ~  Q  (1  -  u»). 

1  1  *0 

In  Figure  3  the  specie  concentrations  are  shown  at  two  downstream  stations 
of  6- 


Next,  the  specie  concentrations  for  the  equilibrium  flow  will  be 
considered.  It  has  already  been  pointed  out  that  over  the  main  part  of 
the  flow  field  the  W  term  is  small  compared  with  the  other  terms  in  the 
specie  equations.  This  statement  immediately  suggests  that  the  specie 
concentration  for  the  equilibrium  flow  may  be  obtained  with  a  perturba¬ 
tion  of  the  pure  diffusion  distribution.  A  general  way  for  doing  this 
would  be  to  first  obtain  C1  with  W  =  0  in  Equation  (20) .  A  better 
approximation  for  Cj  can  then  be  obtained  by  inserting  the  obtained  C^'s 
in  the  production  term,  which  is  essentially  a  product  of  the  various 
C^'8  and  make  a  new  integration.  This  iterative  procedure  can  be  con¬ 
tinued  until  the  C^'s  obtained  are  close  to  those  put  in  the  W  term. 
However,  this  method  is  rather  tedious  and  can  only  be  performed  with  the 
help  of  a  computer.  It  is  far  easier  to  use  the  concept  of  local  chem¬ 
ical  equilibrium  with  the  pure  diffusion  distribution  obtained  from 
Equation  (20)  with  W  =  0.  The  problem  is  then  reduced  to  an  algebraic 
one  and  a  solution  can  br  obtained  based  on  the  thermodynamic  properties 
already  computed  in  the  first  part  of  the  problem.  Let  us  consider  this 
latter  procedure. 


The  solution  of  Equation  (20)  with  W  =  0  and  the  initial  conditions 
Equation  (21)  and  (22)  applied  to  the  initial  equilibrium  composition  is 
first  needed.  Then  at  some  point  in  the  flow  field  we  include  a  sample 
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of  gas  of  the  now  known  diffusion  composition  in  a  small  control  volume 
and  let  this  sample  of  gas  come  to  chemical  equilibrium  at  the  prevailing 
temperature  through  the  following  reactions: 

C02  +  H2  *  CO  +  HaO 

2  lfe  +  02*  2  H20  . 

These  reactions  change  the  diffusion  composition  of  the  gas  in  the  control 
volume.  To  compute  the  new  conposition  we  hnv<^  in  accordance  with  the 
law  of  mass  action,7  for  the  two  reaction:!  above  (p  =  1  atm) 


«n  K  -  (-  £)  +  (-  Jp)  -  (-  lL) 

pl  CO  H20  C02 


(-  Sri  -  *n  r 


CO  *  HgO 


(23) 


H2 


C02  H;, 


in  K  =  2(-  J;)  -  2(-  |{;) 

P2  HaO  m  H2 


"  (_  htI  =  ln 


h2o 


*  r02 


(24) 


Here  K  are  the  reaction  constants,  P  the  free  energy,  and  T.  the 
Pi,2 

molecular  concentration  defined  as 


i 


where  is  the  number  of  gram-molecules  of  specie  i.  The  relation 
between  C^  and  is 


IriMi 


; 


where  Mi  is  the  molecular  weight  of  specie  i. 
obtained  from  standard  gas  tables.8 


The  value  of 


t.L.) 

1  RT  i 


can  be 
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Three  further  equations  for  our  five  unknown  iys  will  be  obtained 
from  continuity  relations.  We  have  first 


Zri  "  1  "  ^non  react  ”  ^Ns*  (2^) 

i 


Furthermore,  in  the  small  control  volume  the  total  amount  of  carbon  to 
oxygen  and  hydrogen  to  oxygen  does  not  change  that  is. 


const . ; 


(26) 


which  may  in  turn  be  translated  to  relations  between  the  concentrations. 
The  five  Equations  (23),  (2k),  (25),  and  (26)  constitute  a  non-linear 
system  for  the  concentrations  at  each  point  which  can  be  solved  by 
numerical  procedures. 

In  Figure  4,  the  equilibrium  specie  concentrations  are  shown  at 
two  downstream  locations,  however,  the  various  curves  are  cut  off  at 
1**1.  It  was  found  in  the  computation  of  the  I^'s  that  for  a  slightly 
larger  value  of  r*  some  of  the  iys  deviated  more  than  50jt  from  their 
diffusion  values.  Therefore,  the  assumption  that  the  W  term  in  Equation 
(20)  is  small  in  this  region  does  not  apply  and  the  result  of  the  simple 
calculation  for  these  r*  values  is  not  very  satisfactory. 


THE  TURBULENT  MIXING  CASE 

As  has  been  shown  by  Pai,5  the  momentum  equations  for  the  turbulent 
and  the  laminar  mixing  processes  can  be  reduced  to  the  same  form,  and 
we  have  the  same  solution  for  the  same  initial  condition's.  We  shall 
briefly  summarize  this  result  and  put  it  in  a  form  suitable  for  the 
present  investigation. 


For  turbulent  Jet  mixing  the  momentum  equation  is1'5: 


on 
?5u  — 
Ox 


(27) 
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where  the  bar  denotes  the  average  value  of  the  variable.  In  general, 
the  eddy  kinematic  viscosity  coefficient  e  may  be  written 

e  =  Cl(|)n  m 

where  n  lies  between  0  and  1,  sx  is  an  empirical  constant  and  L  is  a 
reference  length  which  in  our  case  can  be  chosen  as  rQ.  Both  cx  and  n 
can  be  determined  experimentally  for  the  special  case  under  considera¬ 
tion.  In  terms  of  non-dimensional  quantities  defined  as 


-w  o  u  ~  v  *  ,  o 

p  =  h  u  =  ^  v  =  ^'  €1 


and  the  turbulent  stream  function  defined  such  that 


♦  P  V  .  V  -  **  P  » 

the  momentum  equation  is  reduced  to  (L  =  rQ) 

_ _ — _ x»n  JL(?  ?2r*2u~)  . 
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Here  a  new  independent  variable 


pur 
o  o  o 


y  in  J 

x* _ 

n  + 1 


is  chosen  and  by  letting  cx  Vs  =  1  in  analogy  with  Equation  (11),  Equation 
(Jl)  takes  exactly  the  same  form  as  Equation  (13).  The  resulting  dis¬ 
tribution  for  u  can  be  taken  from  Figures  1  and  2  where  only  £  has  to  be 
replaced  by  =■•  The  scale  factor  between  the  turbulent  and  laminar  dis¬ 
tribution  is  seen  to  be 


Z_  x^‘ 

%  “  n+1  * 


A  similar  argument  can  be  shown  to  hold  also  for  the  specie  concent* ations. 
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DISCUSSION 

The  problem  of  a  multi-component,  hot,  free  Jet  mixing  with  cold 
air  is  very  complex.  To  make  the  problem  tractable  several  simplifying 
assumptions  are  necessary.  With  these  assumptions  in  mind,  however,  it 
cannot  be  expected  that  the  details  of  the  mixing  phenomenon  will  be 
very  well  represented.  Especially  for  the  equilibrium  flow  field  cal¬ 
culation,  for  which  the  initial  temperature  is  rather  high,  the  conjec¬ 
tures  regarding  the  enthalpy  of  the  reacting  mixture  seems  to  be  most 
severe.  However,  by  choosing  a  low  value  of  the  mixture  ratio  the 
chemical  effects  can  be  made  less  important,  so  that  the  simplifications 
employed  can  be  Justified  compared  with  others  inferred  in  the  problem. 

Despite  the  crude  model  it  is  hoped  that  the  main  features  of  the 
mixing  phenomena  are  brought  out. 
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SYMBOLS 


axi-synmetric  coordinates 

stream  function,  defined  in  Equation  (9)  and  (29) 

variable  defined  in  Equation  (12) 

variable  defined  in  Equation  (31) 

component*!  of  velocity  vector 

density 

temperature  °K 

pressure 

Ici<hi  -hi>  -  Ici  t  enthalpy  of  mixture 
i  i 

heat  evolved  in  formation  of  specie  i  at  0°K  per  unit  mass 
free  energy 

mass  concentration  of  specie  i 

molecular  concentration  of  specie  i 

diffusion  coefficient  of  specie  i 

mass  rate  of  production  of  specie  i  per  unit  volume 

viscosity  coefficient 

specific  heat  at  constant  pressure 

heat  conduction  coefficient 

eddy  viscosity  coefficient 


c  H 

Prandtl  number 


cnpDi 

■  ■-  Lewis-Semenov  number 
k 


A,  B 


n 


1 


M 


i 


constants  defined  in  Equation  (7) 
reaction  constant 
constants  defined  in  Equation  (27) 
number  of  gram-molecules 
molecular  weight 


Super-  and  Subscripts 


P‘,  u, . . .  average  value 

P#,  u», . . .  defined  in  Equation  (8) 

P,  u, ...  defined  in  Equation  (28) 

o  refers  to  nozzle  exit 


a 


refers  to  ambient  conditions 
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